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Abstract 

A generic, model-independent method for the analysis of the two-particle short- 
range correlations is presented, that can be utilized to describe e.g. Bose-Einstein 
(HBT or GGLP), statistical, dynamical or other short-range correlation functions. 
The method is based on a data-motivated choice for the zero-th order approxima- 
tion for the shape of the correlation function, and on a systematic determination 
of the correction terms with the help of complete orthonormal sets of functions. 
The Edgeworth expansion is obtained for approximately Gaussian, the Laguerre 
expansion for approximately exponential correlation functions. Multi-dimensional 
expansions are also introduced and discussed. 
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1 Introduction 



Can one model-independently characterize the shape of two-particle correla- 
tion functions? In the present discussion, we address such model-independent 
characterizations of short-range correlations on the level of statistical analy- 
sis. We do not make any theoretical assumptions, neither on the thermal or 
non-thermal nature of the particle emitting source, nor on the negligibility of 
Coulomb and other final state interactions, nor on the presence or the negligi- 
bility of a coherent component in the source, nor on the negligibility of higher 
order quantum statistical symmetrization effects, nor on the negligibility of 
dynamical effects (e.g. fractal structure of gluon-jets) on the short-range part 
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of the correlation function. We simply propose an expansion technique based 
on complete orthonormal sets of functions and an experimental choice for the 
zero-th order approximation to the BECF. 

Although we do not assume here any theoretical knowledge about the short- 
range part of the two-particle correlation function, we assume the following 

experimental properties: 

i) The correlation function tends to a constant for large values of the relative 
momentum Q. 

ii) The correlation function has a non-trivial structure at a certain value of 
its argument. 

The location of the non-trivial structure in the correlation function is assumed 
for simplicity to be close to Q — 0. 

The properties i) and ii) are well satisfied by e.g. the conventionally used two- 
particle Bose-Einstein correlation functions (BECF-s). For an introduction to 
and a review of the recently determined non-ideal features (e.g. non-Gaussian 
shapes in multi-dimensional Bose-Einstein correlation studies), we recommend 
ref. [1]. 

One of the parameters of the fit functions is the intercept parameter A*, that 
carries important information [2,3] on the possible restoration of Ua(X) sym- 
metry in hot and dense matter, when the core-halo picture can be applied [4]. 
However, the intercept parameter A* is known to be particularly sensitive to 
the chosen form of the parameterization of the correlation function. For an ex- 
ponential fitting function, the extrapolated intercept is typically higher than 
that of the Gaussian fit, sometimes even by a factor of two. However, both pa- 
rameters have only a few percent error from the minimization procedure [5,6]. 

It is well known [1], that certain 1 or 2 dimensional correlation functions can 
be fitted with both exponential and Gaussian forms. In case of the NA22 [7] 
and the UAl data [8], the second-order correlation function exhibits a stronger 
than Gaussian or exponential rise at low values of the invariant relative mo- 
mentum Qi. Recently, even singular parameterizations were shown to repre- 
sent well the precision data of the NA22 and the UAl collaborations [7,8]. In 
this formal interpolation of the fitting form to zero relative momen- 

tum would yield infinite value for the intercept. Furthermore, the intercept 
parameter A* of the core-halo picture is rather difficult to measure, as various 
non-ideal effects due to detector resolution, particle mis-identification, reso- 
nance decays, details of the Coulomb and strong final state interactions etc 
may influence this parameter of the fit. One should also mention, that if all of 
these difficulties are properly handled, the intercept parameter A* for like-sign 
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charged bosons is (usually) not larger, than unity as a consequence of quan- 
tum statistics for chaotic sources, even with a possible admixture of a coher- 
ent component. However, attractive final state interactions, fractal branching 
processes of gluon jets, or the appearance of one-mode or two-mode squeezed 
states [9,10] in the particle emitting source might provide arbitrarily large 
values for the intercept parameter. 

Let us emphasize, that our new method is based only on the properties i) and 
ii) of the correlation function. Thus, the method is really model-independent, 
and it can be applied not only to Bose-Einstein correlation functions but to 
any other experimentally determined function, that satisfies properties i) and 
ii). 



2 Shape expansions using complete orthonormal sets of functions 

Let us define the two-particle correlation function as the ratio of the two- 
particle invariant momentum distribution to the product of the single-particle 
invariant momentum distributions [11]: 



where k stands for the three-momentum of the detected particle of mass m 
and energy E = \/m 2 + k 2 . For simplicity, we assume that at large relative 
momenta, the two-particle correlation function tends to 1; i.e. we normalize C 2 
in such a manner that the possibly existing long-range correlations are removed 
by the normalization, which can be achieved by introducing an appropriate 
normalization coefficient into eq. (1). The two-particle correlator is introduced 
as 



We assume that -R 2 (ki, k 2 ) is experimentally determined and it satisfies prop- 
erties i) and ii). Further, we assume that the measured i? 2 (ki,k 2 ) is similar 
to a function w(t), where t denotes certain experimentally defined dimension- 
less combination(s) of relative momentum variables, in 1, 2 or 3 dimensions. 
For example, the experimental data are frequently represented in one dimen- 



sion in terms of t oc Qj = J— (pi — P2) 2 , or in three dimensions in terms of 



t oc (Q , Q z , Q t ) = {E l - E 2 , k M - k 2)Z , yj(k 1<x - k 2 , x ) 2 + (k 1>y - k 2>2/ ) 2 ). We 



define t to be a dimensionless variable, i.e. the relative momentum multiplied 
with certain spatial scale, e.g. t = Q1R1. 






(2) 
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The function w(t) can be considered simultaneously as a zero-th order approx- 
imation to the function R2{t), as well as a certain weight function or measure, 
that is related to an abstract Hilbert space and a complete orthonormal set of 
functions, h n (t). This property of the w(t) function will be the basis of the new 
expansion techique, and the expansion will consist of finding the coefficients 
of the h n (t) functions. For the application of the results to expand the real- 
valued two-particle correlation functions C2(ki,k2), we demand that all the 
h n (t) functions are real- valued functions of the variable(s) t. We discuss those 
situations only, when such a complete orthonormal set of functions exists with 
respect to the measure d/i(t) = dtw(t): 



J dtw(t)h n (t)h m (t) = «J n>m , (3) 

oo 

/W = E/»M*). ( 4 ) 

ra=0 

f n = J dtw(t)f(t)h n (t). (5) 



The Hilbert space H with the measure dtw(t) is spanned by the functions 

{K(t)}™ =0 . 

The important consideration is that we take w(t) to be a function that is simul- 
taneously a measure for an abstract Hilbert space and at the same time, w(t) 
gives a good approximation to the experimentally measured correlator i?2(£), 
so that the first few expansion coefficients would suffice to fully describe the 
measured function. For example, w(t) can be taken as a (multi-dimensional) 
Gaussian function. In case of a one-dimensional Gaussian function, the Her- 
mite polynomials form a complete orthonormal set of functions with respect 
to a Gaussian weight. Then we can obtain a convergent expansion in terms of 
Gaussians and Hermite polynomials as specified in Section 4. Note also that 
the dimensionless variable t can be defined in one, two or three dimensions, 
e.g. t = RiQi or t = (RlQl, RtQt) or t = (RlQl, RsideQside, RoutQout)- The 
multi-dimensional parameterizations will be addressed explicitly in Section 5. 

Let us assume, that the function g(t) = R 2 (t)/w(t) is also an element of the 
Hilbert space H. This is possible, if 



J dtw{t)g 2 (t) = J dt [Rl(t)/w(t)\ < oo, (6) 



i.e. if the experimentally measured correlator, i? 2 (0 approaches the value 
sufficiently fast. With the help of a numerical integration of the experimentally 
determined i?2(0 function and the selected w(t) zero-th order approximation 
to R2(t), the convergence of the integral in the above eq. (6) can be explicitly 
checked. 
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The function g can be expanded as 



g(t) = Y,9nh n (t), (7) 
n=0 

g n = J dtR 2 (t)h n (t). (8) 

Thus the expansion coefficients g n can be determined by a simple numerical 
evaluation of the integral in eq. (8), using the experimentally determined 
values for the correlator R 2 (t). 

From the completeness of the Hilbert space and from the assumption that g(t) 
is in this Hilbert space, i.e. eq. (6) is satisfied, one obtains the following series 
expansion for the correlator R 2 (t): 

oo 

R 2 (t) =w(t)Y / 9nK(t). (9) 

n=0 



It is obvious from this form, that the best choice of the w(t) function is the 
form that is as close to the experimentally determined correlator, as possi- 
ble. For convenience, it is advantageous to choose h (t) = const and use the 
normalization goho(t) = 1, by an appropriate rescaling of the magnitude of 
the function w(t). In such a case, the first few expansion coefficients may be 
sufficient to characterize the correlation function and the correlator R2(t). 

In order to characterize with an independent parameter the strength of the 
two-particle correlator R 2 (t), the expansion-dependent parameter A„, is intro- 
duced, and for convenience, the convention w(t = 0) = 1 is utilized. Thus the 
two-particle BECF can be expanded into the following family of series: 

C 2 (t) = N ll + \ w w(t) jr g„h n (t) \ , (10) 

I n=0 J 

where the core-halo [4] intercept parameter of the correlation function is 

oo 

A* = X w 9nh n (0), (11) 

n=0 



and the coefficients of the expansion, g n can be determined either from nu- 
merical evaluation of eq. (8) or from a fit to the measured data points with a 
particular concrete form of eq. (10). 

The above general example can be utilized to study correlation functions 
in various domains and with various weight functions and related orthogo- 
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nal polynomials, for example the Legrende, associated Legrende polynomials, 
spherical harmonics, Gegenbauer, Chebyshev, Hermite, Laguerre or Jacobi 
polynomials. 

In the next two Sections, we demonstrate the power of the method by working 
out the specific examples for approximately exponential and approximately 
Gaussian correlators, by using w(t) = exp(— t) in the < t < oo domain and 
and w(t) = exp(— 1 2 /2) in the — oo < t < oo domain. As we develop a fitting 
method, we will utilize only the orthogonality of the complete orthonormal set 
of functions {h n (t)}'^ =0 , giving a simple form for H n (t) oc h n (t) and absorbing 
the normalization coefficients of h n {t) into the fit parameters. 



3 Laguerre expansion and exponential shapes 



If in a zeroth-order approximation the correlation function has an exponential 
shape, then it is an efficient method to apply the Laguerre expansion, as a 
special case of eq. (10): 



t — QRl, 
w(t) = exp(— t), 



(JXJ 

J dt exp(-t)L n (t)L m (t)<x5 n>m , 



(12) 
(13) 

(14) 



where the order-n Laguerre polynomial is defined as 

L n (t) = exp(t)— H) n exp(-0. 



(15) 



The general form of eq. (10) takes the particular form of the Laguerre expan- 
sion [6] as: 



a 



(Q)=A/"{ 



l + A L exp(-Qi? L ) 



1 + cMQRl) + ^L 2 (QR L ) + ...] | 



(16) 



The fit parameters are the scale parameters A/", A^, Rl and the expansion 
coefficients Ci, c 2 , ... . The first few Laguerre polynomials are explicitly 
given as 



L (t) = l, 
Li(t)=*-1, 



(17) 
(18) 
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L 2 (t)=t 2 -At + 2, 



(19) 



The Laguerre polynomials are non- vanishing at the origin, hence C 2 {Q = 0) 7^ 
1 + Ax,. The core-halo intercept parameter, A*, is defined by the extrapolated 
value of the two-particle correlation function at Q — 0, see refs. [4,12-14] 
for further details. If the core/halo model is applicable, and effects of partial 
coherence in the particle emitting source can be neglected, the core-halo in- 
tercept parameter A* measures the squared fraction of bosons emitted directly 
from the core, an important physical observable that is related to the degree of 
partial restoration of £7^(1) symmetry in hot and dense hadronic matter [2,3], 
as a partial £7^(1) symmetry restoration leads to enhanced production of rj' 
mesons and suppression of the core fraction of directly produced pions at low 
values of transverse momentum. 

The physically significant core-halo intercept parameter A*, and its error <5A* 
can be obtained from the parameter Ax of the Laguerre expansion as 



A* = Ax[l - ci + c 2 - ...], 
5 2 K = 5 2 \ L \l + c\ + c\ + ...1 + \\ \5 2 d + 5 2 c 2 + ... 



(20) 
(21) 



Any Bose-Einstein correlation function can be expanded into a convergent 
Laguerre expansion of the form of eq. (16), provided that 



00 

J dt R 2 2 {t) exp(-K) < 00, 



(22) 



where t = QRl stands for the dimensionless scale variable. In principle, the 
left-hand side of this unequality can be evaluated directly from the experi- 
mental data. 



4 Edgeworth expansion and Gaussian shapes 



If, in a zeroth-order approximation, the correlation function has a Gaussian 
shape, then it is an efficient method to apply the Edgeworth expansion to 
characterize the deviation from the Gaussian form in the following manner: 



* = V2QR E , 
w(t)=exp(-t 2 /2), 



J dtexp(-t 2 /2)H n (t)H m (t)^5 n , m , 



(23) 
(24) 

(25) 
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where the order-n Hermite polynomial is defined as 

H n (t) = exp(t 2 /2) exp(-t 2 /2). 



(26) 



The general form of eq. (10) takes the particular form of the Edgeworth ex- 
pansion [15,5,6] as: 



C 2 (Q)=N{l + X E exp{-Q 2 R 2 E ) x 



3! 



1 ~H 3 (V2QR E ) + ^H 4 (V2QR E ) + ... 



(27) 



The fit parameters are the scale parameters Af, X E , Re and the expansion 
coefficients k 3 , k 4 , ... , that coincide with the cumulants of rank 3, 4, of 
the correlation function. The first few Hermite polynomials are listed as 



i?i (*)=*, 
H 2 (t) = t 2 -l, 
H 3 (t) = t 3 -3t, 
H±(t) = t 4 -6t 2 + 3, 



(28) 
(29) 
(30) 
(31) 



The physically significant core-halo intercept parameter A* can be obtained 
from the Edgeworth fit of eq. (27) as 



A* — X E 

f2 x _ X 2 



/t 4 



S 2 \, = S 2 X E + (k 4 SX e + X e Sk 4 Y/64 + ... . 



(32) 
(33) 



Any Bose-Einstein correlation function can be expanded into a convergent 
Edgeworth expansion, if 



J dt Rj(t) exp(+t 2 1 2) < oo. 



(34) 



This latter requirement can be checked experimentally, if necessary. This ex- 
pansion technique was applied in the conference contributions [5,6] to the 
AFS minimum bias and 2-jet events to characterize successfully the deviation 
of data from a Gaussian shape. It was also successfully applied to characterize 
the non-Gaussian nature of the correlation function in two-dimensions in case 
of the preliminary E802 data in ref. [6], and it was recently applied to char- 
acterize the non-Gaussian nature of the three-dimensional two-pion BECF in 
e + + e~ reactions at LEP - 1 [16]. 
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5 Multi-dimensional expansions 

Genuine mult i- dimensional expansions can be obtained by using multi-dimensional 
complete orthogonal sets of functions. The simplest two-dimensional example 
is the study of angular correlations on the surface of the unit sphere (e.g. to 
study correlations of jets within jets in QCD). In this case, t = (0,<f>), the 
measure is w(t) = sin(#) and the domain is (0 < 6 < ir, < (f) < 2ir). The 
orthogonal polynomials are the well-known spherical harmonics Y l m (9, 0) that 
satisfy 

it 2tt 

J d6 sin(0) J d<t> ¥1(9, (j>)YT(e, 0) oc 6 k j6 m , n (35) 



and the angular correlations can be expanded in a series as 

C 2 (9,</>) =l + A y sin(0)[l + f; £ c lm Yr(e,<f>)}. (36) 

1=1 m=—l 

This example should be sufficient to demonstrate how the method can be gen- 
eralized to higher dimensions. A systematic treatise of the multi-dimensional 
realizations of the general expansion method as described in Section 2 is be- 
yond the scope of the present paper. 

In the following we consider the special case of factorized multi-dimensional 
distributions. A more profound technique can be obtained based on the special- 
ization of the general arguments of section 2 to the case of multi-dimensional 
BECF-s. The two-dimensional Edgeworth expansion and the interpretation 
of its parameters is described in the handbook on mathematical statistics by 
Kendall and Stuart [17]. 

If the two-particle BECF can be factorized as a product of (two or more) func- 
tions of one variable each, then e.g. a Laguerre or an Edgeworth expansion 
can be applied to the multiplicative factors - functions of one variable, sepa- 
rately. This method was applied recently to study the non-Gaussian features 
of multi-dimensional Bose- Einstein correlation functions in refs. [5,16]. 

Let us restrict ourselves to the situation when the mult i- dimensional Fourier- 
transformed emission function factorizes in certain set of variables. For ex- 
ample, such a boost-invariant form of factorization is assumed in the Buda- 
Lund parameterization of the two-particle Bose-Einstein correlation functions 
in refs. [18,19]. 
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A simple expression for the two-particle correlation functions is obtained for 
longitudinally approximately boost-invariant, cylindrically symmetric particle 
emitting sources, in the Buda-Lund picture, if we assume that the emis- 
sion function factorizes as a product of an effective proper-time distribution, 
a space-time rapidity distribution and a transverse coordinate distribution 
where 77 = 0.5 \og[(t + r z )/(t — r z )\ is the space-time rapidity, r = ^t 2 — r\ is 
the longitudinally boost-invariant proper-time distribution [18,19]. 

With the help of a small source size (or large relative momentum) expansion, 
the two-particle Bose-Einstein correlation function can be written into the 
following boost-invariant, factorized Buda-Lund form: 



C (Ak K\ - 1 + A (K) |AW - >|2 l 6 -Wll)l 2 feglll! r 3 7) 
C 2 (Ak,K)-l + X.(K) |#t(())|2 |&(Q)|2 m? . (37) 

Here, the Fourier-transformed proper-time , space-time rapidity and transver- 
sal coordinate distributions, H*(Q=), G*(Q\\) and I*(Q±) can be of power-law, 
exponential, Gaussian or other types, corresponding to the underlying space- 
time structure of the particle emitting source, hence it is natural to apply 
the one-dimensional Edgeworth or the Laguerre expansions to any of these 
factors. 

The invariant temporal, parallel and perpendicular relative momentum com- 
ponents are defined, respectively, as 



Q = = Qocosh[r]] - Q z smh[i]} = Q ■ n, (38) 
Q» = Q sinh[?j] — Q z coshfrj] = Q x n, (39) 



Q± = ^Ql + Q 2 y , (40) 

Q 2 = Q-Q = (Q = ) 2 - (Qn) 2 - Q 2 ± . (41) 

In the above, a ■ b = a M 6 M = ao&o — ab = aob — a x b x — a y b y — a z b z stands 
for the inner product of four-vectors. The direction of the normal-vector n is 
characterized as n M = (cosh [77], 0, 0, sinhfrj]), by a mean space-time rapidity 
r\ [18,19] in the LCMS frame [20]. Hence this direction is boost-invariant and 
fj is one of the fit parameters. As n ■ n = +1, this normal- vector points into 
a time-like direction and Q -n = Q = is an invariant time-like component of 
the relative momentum. The relative momentum component Q\\ is parallel to 
the longitudinal direction (which is the thrust axis in jet physics, or the beam 
direction in heavy ion physics), and Qn is invariant to boosts along this direc- 
tion. Finally, Qj_ is the remaining perpendicular (or transverse ) component 
of the relative momentum, it is also invariant for longitudinal boosts. 
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In eq. (37) the Fourier-transformed distributions are defined and explained 
in greater detail in refs. [18,19]. Although these distributions can be theo- 
retically evaluated from an assumed shape of the emission function S(x,k), 
for example see ref [18,19], here we are interested in the model-independent 
characterization of the Bose-Einstein correlation functions, exclusively. 

From a three-dimensional analysis of two-pion Bose-Einstein correlation data 
at LEP1 [16] we know that the BECF is approximately a factorized Gaus- 
sian. Thus, a good candidate to characterize the non-Gaussian nature of these 
correlations in 3 dimensions, in a boost-invariant manner, is the following 
Edgeworth expansion of the factorized form of the Buda-Lund correlation 
function: 



C 2 (Ak, K) = l + \ E exp(-Q 2 = Ri - Qjjfljj - Q{Ri) x 

1 + ^H 3 (V2Q = R = ) + ^H,(V2Q = R = ) + ... 



x 



1 + 
1 + 



H 3 (V2Q\\R\\) + ^H A (V2Q\\R\ { 



+ ... 



x 



«3,|| 

3! " ov v — *ir-ii' 1 4! 
^H 3 (V2Q ± R ± ) + ^H 4 (V2Q ± R ± ) + ...] . (42) 



What are the fit parameters in eq. (42)? We have 5 free scale parameters 
for cylindrically symmetric, longitudinally expanding sources, and three se- 
ries of shape parameters. The scale parameters are \e, R=, R\\, R± and rj, 
that characterize the effective source at a given mean momentum, by giv- 
ing the vertical scale of the correlations, the invariant temporal, longitudinal 
and transverse extensions of the source and its invariant direction, which is 
the space-time rapidity of the effective source in the LCMS frame (the frame 
where ki jZ + k 2 , z = 0, [20]). The three series of shape parameters are /c 3 _, 

«4,= , ••• , «3,||, «4,||, ••• , «3,±, «4,±, ••• • 

In principle, each of these parameters may depend on the mean momentum K. 
At any fixed value of the mean momentum K, the above free parameters of the 
invariant Buda-Lund correlation function can be fitted to data, without any 
a priori assumption on the shape of the source emission function, based only 
on the cylindrical symmetry and the boost-invariant factorization property of 
the particle emitting source. 



6 Application of the results 



As the Laguerre and the Edgeworth expansions of Bose-Einstein correlation 
functions were introduced with quite some success in refs. [5,6] to characterize 
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the AFS 2-jet and the AFS minimum-bias data in 1 dimension, to characterize 
the non-Gaussian features of the E802 preliminary Si + Au — > 7r + + 7r + + X 
data in two dimensions [6] , and they were reported to successfully characterize 
the non-Gaussian features of the L3 two-pion BECF at LEP1 in two and three 
dimensions, [16] (although using a non-invariant formulation), we think that 
there is no need to further demonstrate the applicability of the Edgeworth 
expansion and its multi-dimensional generalizations in correlation studies in 
high energy physics. 

However, the applications of the Laguerre expansions were not yet presented 
in the literature as far as we know. Hence, let us demonstrate in Figure 1 
the power of the Laguerre expansions to characterize two well-known, non- 
Gaussian correlation functions: the CiiQi) correlation functions as determined 
by the UAl and the NA22 experiments [7,8]. Note that the invariant momen- 
tum difference Q 2 is binned in logarithmic scale, and the resolution of both 
experiments goes down to about 30 MeV, a very small scale as compared to 
the typical 200 MeV half-width of the BECF in particle physics. 

We note, that Fig. I. is just an illustration of the power of Laguerre expansions, 
however, the authors do not intend to replace the power-law description of 
these correlation functions by the Laguerre expansion, and the authors are 
aware of the fact that a power-law fit gives a similarly good x 2 /NDF with 
smaller number of fit parameters. 

Table 1 



Best fits to UAl and NA22 two-particle correlations using a Laguerre expansion 





UAl 


NA22 


Parameter 


Value 


Error 


Value 


Error 


N 


1.355 


± 0.003 


0.95 


± 0.01 


Xl 


1.23 


± 0.07 


1.37 


± 0.10 


Rl [fin] 


2.44 


± 0.12 


1.35 


± 0.14 


Cl 


0.52 


± 0.03 


0.63 


± 0.06 


C2 


0.45 


± 0.04 


0.44 


± 0.06 


X 2 /NDF 


41.2/41 


= 1.01 


20.0/34 = 


0.59 



From Table 1, one can determine that the quality of the Laguerre expan- 
sion fits is statistically acceptable both in case of the UAl and in case of the 
NA22 data. We can also determine the core-halo model intercept parameter, 
A, = 1.14 ± 0.10 (UAl) and A* = 1.11 ± 0.17 (NA22). As both of these values 
are within errors equal to unity, the maximum of the possible value of the 
intercept parameter A* in a fully chaotic source, we conclude that either there 
are other than Bose-Einstein short-range correlations observed by both collab- 
orations, or in case of this measurement the full halo of long-lived resonances is 
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resolved [4,13]. This is very interesting and the distribution of long-lived reso- 
nance fractions and decay times may be related to the approximate power-law 
shape of the fitted correlation functions [22]. The full resolution of the halo 
is theoretically possible only because the rj'-s and 77-s were either produced in 
a negligible number, or because their effects were corrected for in the Monte- 
Carlo evaluation of the background distribution. 



7 Summary 



In this paper we described a general method to characterize, in a model and 
theory independent manner the shape of the two-particle Bose-Einstein cor- 
relation functions with the help of mutually orthogonal sets of functions, as- 
suming that the measure of the abstract Hilbert space can be identified with 
the approximate shape of the correlation function. For one- dimensional, ap- 
proximately Gaussian correlation data, the Edgeworth expansion technique is 
obtained. For one- dimensional, approximately exponential shapes, we obtained 
the expansion in terms of the Laguerre polynomials. We generalized these ex- 
pansion techniques to two and three dimensions. We also demonstrated that 
the method works well: we characterized the deviation of the two-particle 
Bose-Einstein correlation function of the UAl and the NA22 collaborations 
from the exponential shape with the first two Laguerre coefficients. 
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Laguerre expansion fit 




Fig. 1. The figures show Z?| which is proportional to the two-particle Bose-Einstein 
correlation function, as measured by the UA1 and the NA22 collaborations. The 
dashed lines stand for the exponential fit, which clearly underestimates the measured 
points at low value of the squared invariant momentum difference (note the 
logarithmic horizontal scale). The solid lines stand for the fits with the Laguerre 
expansion method, which is able to reproduce the data with a statistically acceptable 
X 2 /NDF. The fit results are summarized in Table 1. 
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